Isolated Dirac cone on the surface of ternary tetradymite-like topological insulators 
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We have extended our search for topological insulators to the ternary tetradymite-like compounds 
M2X2Y (M=Bi and Sb; X and Y = S, Se, and Te), which are variations of the well-known binary 
compounds Bi2Se3 and Bi2Te3- Our first-principles computations suggest that five existing com- 
pounds are strong topological insulators with a single Dirac cone on the surface. In particular, 
stoichiometric Bi2Se2S, Sb2Te2Se, and Sb2Te2S are predicted to have an isolated Dirac cone on 
their naturally cleaved surface. This finding paves the way for realizing the topological transport 
regime. 



Topological insulator is a very distinct state of quan- 
tum matter with a bulk-insulating gap of spin-orbit cou- 
pling (SOC) originPE] The Dirac fermions on the sur- 
face of a topological insulator possess topological quan- 
tum Berry phase and a definite spin chirality, which pro- 
tects them from back-scattering and localization, result- 
ing in counterpropagation of opposite spin states without 
dissipation. These are the key electronic ingredients for 
applications to topological quantum computing } 12 * 13 * and 
low-power spintronics^ devices for which, however, the 
dissipationless surface states must be in the topological 
transport regime, i.e. to have an isolated Dirac cone fully 
separated from the bulk bands and the Fermi level to lie 
at the Dirac point. 15 The existing topological insulating 
materials fail in this regard due to their material-specific 
complications such as interference from numerous sur- 
face states^, shieldi ng of the Dirac cone either inside 
the bulk- valence bancP^U^ or the conduction bands^, or 
the inv olvem ent of crystal distortion in fabricating the 
materiaPHU. 

Recently, it was shown experimentally that an isolated 
Dirac cone can be achieved by tuning the Fermi level 
with appropriate hole doping, but the non-stoichiometric 
crystal structure and the requirement of surface deposi- 
tion generally makes this approach unsuitable for most 
practical applications lillisj Q n ^ ne other hand, solid solu- 
tions in the tetradymite-like compounds provide a plat- 
form for engineering topological surface states. For exam- 
ple, theoretical work of Zhang et al. 22 indicates the pres- 
ence of a topological phase transition in Sb2(Tei_ :c Se ;E )3 
for < x < 1. Interestingly, since Se atoms preferen- 
tially occupy the central site in the crystal, at #=1/3 
the stoichiometric compound Sb2Te2Se can admit an or- 
dered phase. Our recent work on topological insula- 
tor Bi2Te2Se reveals that the linewidth of the topolog- 
ical surface states in angle-resolved photoemission spec- 
troscopy (ARPES) is narrow, indicating that disorder ef- 
fects are suppressed in these compounds.^ Ren et alP^ 
have shown that Bi2Te2Se is the best material to date for 



studying the topological surface states due to its highest 
bulk resistivity among all known topological insulators. 

Here, we report first-principles computations on five 
existing tetradymite-like compounds Bi2Te2Se, Bi2Te2S, 
Bi2Se2S, Sb2Te2Se and Sb2Te2S to show that these ma- 
terials host salient topological insulating features in their 
stoichiometric crystal structures. Our important findings 
include that the last three of the five aforementioned com- 
pounds harbor an isolated Dirac cone with Dirac points 
which lie fully within the gap, allowing these states to 
reach the long-sought territory of the topological spin- 
transport regime^ 5 -^ where the in-plane carrier trans- 
port would have a purely quantum topological origin. 

Topological insulators in two- and three-dimensions 
have been predicted theoretically and observed experi- 
mentally as well to display quantum spin Hall effecP™ 
and strong topological behavior due to the combined 
effects of relativistic Dirac fermions and quantum- 
entanglement^"^. For systems with inversion symme- 
try, a band inversion occurs due to a finite value of the 
spin-orbit coupling. This results in a bulk-insulating 
gap which ensures the existence of metallic surface states 
within the gap. 16 In bismuth or antimony metals, numer- 
ous surface states interfere with each other and an odd- 
number of Dirac cones can be achieved by making com- 
posite alloys such as Bii^Sb^P^In other novel classes of 
topological insulators, the Dirac node appears either be- 
low the valence band maximum as in Bi2Te3 or Sb2TesP^ 
or above the conduction band minimum as in TlSbTe^Pl 
Topologically nontrivial materials such as half-Heuslei^ 
and Li2AgSb 21 series are semi-metals, which require lat- 
tice distortion to be tuned into the topological insulat- 
ing state. For these reasons, it is difficult to electrically 
gate these materials for the manipulation and control of 
charge carriers for realizing a device. Therefore, a viable 
route to search for the isolated Dirac node in topological 
insulating materials is to find a chemically stable stoichio- 
metric composition in an appropriate class of materials. 

With this motivation, we consider in this article five ex- 
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FIG. 1. Tetradymite crystal structure and topological 
transport regime, a, The crystal structure of tetradymite 
M2X2Y. M, X, and Y atoms are denoted by red, blue, and 
grey balls, respectively. Diagrams in b and c illustrate surface 
band structures near the T-point for stoichiometric Bi2Te3 
and Bi 2 Se2S. 



isting stoichiometric tetradymite-like ternary compounds 
in the structure M 2 X 2 Y , namely, Bi 2 Te 2 Se^, Bi 2 Te 2 ^, 
Bi 2 Se 2 Sp§] Sb 2 Te 2 SeP3 and Sb 2 Te 2 ^D. Tetradymite 
structure is with a rhombohedral unit cell of the space 
group R3m. The commonly invoked hexagonal cell 
shown in Fig. la consists of three quintuple layers. The 
stacking order of two consecutive units may be repre- 
sented as -X-M-Y -M-X--X-M-Y -M-X-. 
The layer of Y atoms is in the center of the unit. The 
natural surface termination is between the two X-atom 
layers. The known topological insulators Bi 2 Se3, Bi 2 Te3, 
and Sb 2 Te3 have the same crystal structure with X = Y. 

Fig. lb shows a schematic illustration of the bulk and 
surface band structures of Bi 2 Se3, Bi 2 Te3 and Sb 2 Te3 in 
which the surface Dirac point lies below the valence band 
maximum (VBM). In our earlier experimental work, we 
lowered the Fermi level (Ef) of Bi 2 Se3 into the bulk 
bandgap by substituting trace amount of Ca 2+ for Bi 3+ 
with Ca acting as a hole doner to the bulk states. In or- 
der to lift the Dirac point above the VBM, we deposited 
N0 2 on the surface. An isolated Dirac cone is then ob- 
tained as illustrated in Fig. lc. In the present study, 
we predict new topological insulators Bi 2 Te 2 Se, Bi 2 Te 2 S, 
Bi 2 Se 2 S, Sb 2 Te 2 Se and Sb 2 Te 2 S, among which the latter 
three have a Dirac point located in the bulk gap above 
the VBM as shown in Fig. lc. Therefore, these three 
compounds are well-suited as materials with an isolated 
Dirac cone without requiring tuning via surface deposi- 
tion. With proper electron/hole doping control or gat- 
ing, one can position the Ef at the Dirac point and the 
material can be in the transport regime for potential de- 
vice applications. We did not attempt to tune the elec- 
tronic structure of the studied compounds us ing vi rtual 
crystal 3 -^ 3 - 3 -! or other first-principles approache d 34 * 35 *. 

The first-principles band structures were computed on 
the basis of experimental lattice datsP^H^H within the 
framework of the density-functional theory (DFT) us- 
ing the generalized-gradient approximation (GGA). Our 
results, shown in Fig. 2, predict that all these materials 
naturally host a topological insulating phase with band 
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FIG. 2. Bulk electronic structures and band- 

inversion. Bulk electronic structures of Bi 2 Te2Se (a), 
Bi 2 Te 2 S (b), Bi 2 Se 2 S (c), Sb 2 Te 2 Se (d) and Sb 2 Te 2 S (e). 
The grey shaded areas in panels a-e highlight the inverted 
band structure as illustrated in panel f. 



gaps in the range of 224 meV to 297 meV. A band inver- 
sion occurs at the T-point as highlighted with grey shad- 
ings. As illustrated in Fig. 2f, the conduction band (red 
solid line) and the valence band (blue dashed line) cross 
each other near the T-point, so that at this point the oc- 
cupied state possess the character of the conduction band 
and the unoccupied state that of the valence band. Since 
the crystal has inversion symmetry, we have applied band 
parity analysis 5 to evaluate the value of the topological 
invariant Z 2 . All five compounds are thus found to be 
strong topological insulators with Z 2 = —1. This comes 
about due to band inversion at the T-point through which 
the Z 2 topological invariant picks up an extra factor of 
-1 compared to bands without SOC. In short, the pres- 
ence of an energy gap throughout the Brillouin zone with 
a band inversion only at the T-point makes this class of 
materials topological insulators. 

The calculated DFT-GGA surface bands are shown in 
Fig. 3. A pair of surface bands with opposite spin di- 
rections are degenerate at the zone center and form the 
Dirac-cone-like dispersion. All five compounds exhibit a 
single Dirac cone inside the bulk band gap. For Bi 2 Te 2 Se 
and Bi 2 Te 2 S, the Dirac point lies below the VBM. This is 
the same as the case of stoichiometric binary Bi 2 Te3 com- 
pound in which bulk scattering channels at the energy of 
the Dirac point prevent interesting quantum phenomena 
in the topological transport regime to be realized. The 
most interesting results are found in Bi 2 Se 2 S, Sb 2 Te 2 Se, 
and Sb 2 Te 2 S where an isolated Dirac cone is obtained. 
Here, the Dirac point lies above the VBM and below the 
conduction bands. When Ef is tuned to the Dirac point, 
the density of states approaches zero and there are no 
other scattering channels. Therefore, Bi 2 Se 2 S, Sb 2 Te 2 Se, 
and Sb 2 Te 2 S would be stoichiometric compounds which 
are wll-suited for investigating the topological transport 
regime. 
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FIG. 3. Surface states and spin-texture of the Dirac 
cone. Surface band structures of Bi 2 Te2Se (a), Bi 2 Te2S (b), 
Bi 2 Se 2 S (c), Sb 2 Te 2 Se (d) and Sb 2 Te 2 S (e). Red lines are 
the surface states and the blue shaded areas are projected 
bulk states. Surface momentum is defined in panel f showing 
the hexagonal surface Brillouin zone. Green arrows in panel f 
denote the spin direction of the upper Dirac cone at constant 
energy (blue circle). 



In conclusion, we predict that stoichiometric 
tetradymite-like compounds Bi 2 Te2Se, Bi 2 Te2S, Bi 2 Se2S, 
Sb2Te2Se, and Sb2Te2S are strong topological insulators 
with large enough band gap for room temperature 
applications. The isolated Dirac cone on the surface of 
Bi 2 Se2S, Sb2Te2Se, and Sb2Te2S will be particularly 
important for studying topological transport properties 
on a bulk material without the necessity of surface 
deposition.^ Experiments with magnetic scattering 
techniques should be able to uniquely access the helical 
spin-texture of the isolated surface states which are 



essential in device applications. The stoichiometric 
crystals of the investigated compounds are very easy 
to produce and manipulate for manufacturing de vices 
for microchips spintronics^ and informatio n^ 2 * 13 * 
technologies. Furthermore, the full exposure of the topo- 
logical transport regime for dissipationless spin-current 
with the unique advantage of tunable surface states 
is suitable for studying novel topological phe nomenon 
such as intrinsic quantum spin Hall effect pEEU universal 
topological magneto-electric effect 37 and anomalous 
half-integer quantization of Hall conductance^ 3 -^, and 
the quantum electrodynamical phenomena such as 
the image magnetic monopole induced by an electric 
charge 39 or Majorana fermions induced by the proximity 
effects from a superconductor 12 . 
Methods: 

First-principles band calculations were performed with 
the linear augmented-plane-wave (LAPW) method us- 
ing the WIEN2K package 40 within the framework of the 
density functional theory (DFT). The generalized gra- 
dient approximation (GGA) was used to describe the 
exchange-correlation potential. 41 Spin orbital coupling 
(SOC) was included as a second variational step using 
a basis of scalar relativistic eigenfunctions. The sur- 
face was simulated by placing a slab of 6 quintuple layes 
for Bi2Te2Se and Bi2Te2S. For Bi2Se2S, Sb2Te2Se and 
Sb2Te2S we used 12 quintuple layers to ensure conver- 
gence with the number of layers. 
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